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Abstract 

The use of partial geometries to construct parity-check matrices for LDPC codes has resulted in 
the design of successful codes with a probability of error close to the Shannon capacity at bit error 
rates down to 10 -15 . Such considerations have motivated this further investigation. A new and simple 
construction of a type of partial geometries with quasi-cyclic structure is given and their properties 
are investigated. The trapping sets of the partial geometry codes were considered previously using the 
geometric aspects of the underlying structure to derive information on the size of allowable trapping sets. 

This topic is further considered here. Finally, there is a natural relationship between partial geometries 
and strongly regular graphs. The eigenvalues of the adjacency matrices of such graphs are well known 
and it is of interest to determine if any of the Tanner graphs derived from the partial geometries are 
good expanders for certain parameter sets, since it can be argued that codes with good geometric and 
expansion properties might perform well under message-passing decoding. 

Index Terms 

LDPC codes, partial geometries, strongly regular graphs, protographs, expander graphs 

I. Introduction 

Partial geometries play an important role in the construction of low-density parity-check 
(LDPC) codes lfl6l . |[32l . l(33l which currently give the most promising coding technique for 
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error control in communication and data storage systems due to their capacity-approaching per¬ 
formances and practically implementable decoding algorithms. Partial geometries are members 
of a broad class of combinatorial configurations with geometric properties referred to as finite 
geometries. The first classes of LDPC codes based on partial geometries were constructed based 
on Euclidean and Projective geometries over finite fields ll24ll . These classes of finite geometry 
LDPC codes have an abundance of algebraic and geometric structures and perform well with 
iterative decoding algorithms lf24l . [f32ll . lf33ll . The construction of LDPC codes given in lf24l 
was later generalized in different directions ll22l . lt42ll . Il43l . Il53l . which resulted in several large 
classes of finite geometry LDPC codes. Codes based on the more general partial geometries 
were presented in ll20l OTI . Il28l . j48l . In a recent paper lITTIl . it was shown that diverse classes 
of LDPC codes that appear in the literature are actually partial geometry codes although their 
construction methods were not based on geometric notions. 

In this paper, aspects of partial geometries and their use in coding theory are considered. 
A characterization of a special category of partial geometries realized from an array of cyclic 
permutation matrices is given. The graph representation of a partial geometry is also described. 
Two new classes of partial geometries are constructed, one from prime fields and the other one 
from cyclic subgroups of prime orders of finite fields. New classes of quasi-cyclic (QC) LDPC 
codes are constructed based on these new constructions of partial geometries. The problem of 
determining the sizes of trapping sets of such geometric codes, initiated in fTO . is continued 
here for both the partial geometries and the subclass of generalized quadrangles (GQs). Linally, 
since the eigenvalues of the adjacency matrices of all the geometric objects under investigation 
are known, the expansion properties of them can be determined and are summarized to consider 
the possibility that they may play a role in suggesting promising candidates for LDPC codes for 
further investigation. Such candidates would have to be verified by simulation. 

II. Definitions, Concepts and Structural Properties of Partial Geometries 

Consider a system composed of a set N of n points and a set M of m lines where each line 
is a set of points. If a line L contains a point v, we say that v is on L and that L passes through 
v. If two points are on a line, then we say that the two points are adjacent and if two lines pass 
through the same point, then we say that the two lines intersect, otherwise they are parallel. 
The system composed of the sets N and M is a partial geometry 0, 0, [HBJ if the following 
conditions are satisfied for some fixed integers 7 > 2, p > 2 and 5 > 1 : 
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1) Any two points are on at most one line; 

2) Each point is on 7 lines; 

3) Each line passes through p points; 

4) If a point v is not on a line L, then there are exactly <5 lines, each passing through v and 
a point on L. 

Such a partial geometry will be denoted by PaG( 7 , p. 5) and 7 , p and 5 are called the parameters 
of the partial geometry. The parameter 5 is called the connection number of the geometry. 

A simple counting argument [[51 shows that the partial geometry PaG( 7 , p, 5) has exactly 

ri = p({p~ l ){ 7-!) + <*)/£ (!) 

points and 

m = 7((p - 1)(7 - !) + 8)/6 ( 2 ) 


lines. 

If v and v' are adjacent points, then there are exactly 7 A -E/O — 7 — <5 — 1 points, such that 
each of these points is adjacent to both v and v'. On the other hand, if v and v' are not adjacent, 
then there are exactly 7 ^ points, such that each of these points is adjacent to both v and v'. 
Each point v is adjacent to (5 + l)(p — 1) other points. Three adjacent non-colinear points form 
a triangle. It follows from the above adjacency property that two adjacent points v and v' are 
on 75 + p — 7 — 5 — 1 triangles. As will be discussed later, such a triangle in the (point-point) 
adjacency matrix of the graph of the (point-point) partial geometry leads to a cycle of length 6 
in the corresponding Tanner graph ll44l . assuming 5 > 1. 

Well known examples of partial geometries are Euclidean and projective geometries [0, li34l 
over finite fields. If 5 = 7 — 1, the partial geometry PaG( 7 , p, 7 — 1) is called a net which consists 
of n = p 2 points and m = 7 p lines. Each point v not on a line L is on a unique line which is 
parallel to L. Equivalently, there is a unique line IJ that intersects v and L. The set of m = 7 p 
lines in the net PaG( 7 , p, 7 — 1) can be partitioned into 7 classes, each consisting of p lines, 
such that all the lines in each class are parallel, any two lines in two different classes intersect, 
and each of the n = p 2 points is on a unique line in each class. These classes of lines are called 
parallel bundles. A two-dimensional Euclidean geometry (or affine geometry) Il34il is a net. 

For every point v in PaG( 7 ,p, 7 — 1), there are exactly 7 lines that intersect at v, i.e., lines 
that pass through v. These lines are said to form an intersecting bundle at v, denoted by A(v). 
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Notice that v is on every line in A(u), there are exactly 7 (p — 1 ) points, each is on a unique 
line in A(t>), and all the other n — y(p — 1) — 1 points in PaG( 7 , p, 7 — 1) are not on any line 
in A(v). If 8 = p, then every point in PaGlv, p, p) is adjacent to v since every point is on a 
line in A(n). In this case, any two points in PaGfy, p, p) are connected by a line. Examples for 
which 5 = p include two-dimensional Euclidean and projective geometries (also called affine 
and projective planes) lf34l . 

Denote the points and lines in PaG( 7 , p, 5) by v 0 , Vi ,..., v n -\ and L 0 , Li,..., L m _ 1 , respec¬ 
tively. Then, N = {n 0 , 17 ,..., n n _i} and M = {L 0 , L\,, L m _ ,}. Algebraically, a partial 
geometry PaG( 7 , p, 5) with n points and m lines is commonly represented by a m x n matrix 
Hp aG = [hi,j]o<icm, o <j<n with 0 and 1 entries whose rows and columns correspond to the lines 
and points of PaG( 7 ,p, 8), respectively. The rows are labeled from 0 to m — 1 (or by the lines 
{L 0 , L i,..., L m _|} in this order) and the columns are labeled from 0 to n — 1 (or the points 
{t>o,ui, • • • ,v n -i} in this order). The entry h,, 3 is 1 (i.e., h t j = 1) if and only if the point re¬ 
labeled by j (called the j-th point) is on the line L, labeled by i (called the z-th line); otherwise, 
hi,j = 0. This matrix H Pa o is called the line-point adjacency matrix of the partial geometry 
PaG( 7 , p, 5) which shows the incidence relationship of the lines and points of PaG( 7 , p, 5). The 
i-th row of H PaG is called the incidence vector of the i-th line L, . The transpose Hp aG of H PaG is 
also an adjacency matrix of PaG( 7 , p, 8) called the point-line adjacency matrix of PaG( 7 , p, 8). 

Graphically, a partial geometry PaG( 7 , p, 5) can be displayed by an adjacency graph Giy, p, 5) 
with its n points represented by n nodes or vertices, labeled by v 0 ,Vi,... ,n n -i- Two vertices 
Vj and Vk are connected by an edge (ly. zy) if and only if they are on the same line. A line 
L = {vj 0 , v 3l ,..., Vj p _ } ] in PaG( 7 , p, 8) is represented by a sequence of p — 1 connected edges, 
( y j 0 , Cjj), (vji, v j 2 ),..., (vj p _ 2 , Each node v jk in G(y,p,8) has 7 edges incident with it 

and 7 is called the degree of v Jk . It is the number of lines in PaG( 7 , p. d) that intersect at the 
point v jk . 

The partial geometry PaG( 7 ,p, 5) can also be represented by a bipartite graph PTll . denoted 
by G(V, C), which is more commonly used in coding theory. It consists of two sets of nodes, 
denoted by V ( the set of variable nodes (VNs)) and C (the set of check nodes (CNs)), of size 
n and m, respectively. The n nodes in V, denoted by v 0 , 17 ,..., v n _ 1 , represent the n points in 
PaG( 7 , p, <5) with v 3 representing the point v 3 and the m nodes in C, denoted by c 0 , c 1: ..., c m _ 1 , 
represent the m lines in PaG( 7 , p, 5) with c, representing the line L,. A node v 3 in V is connected 
to a node c, in C if and only if the point v 3 in PaG( 7 , p, 5) represented by the node v 3 is on 
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the line L , in PaG( 7 , p. 5 ) represented by the node c % . This bipartite graph Q(Y. C) is called the 
Tanner graph ll44l associated with the line-point adjacency matrix H PaG of the partial geometry 
PaG( 7 ,p, <5), commonly used in the study of LDPC codes in coding theory [|44|. We call this 
bipartite graph a ( 7 , p )- biregular bipartite graph since the VNs on the left are regular (have the 
same degree 7 ) and the CNs on the right are regular (of the same degree p ). The associated (or 
the adjacency) matrix H PaG is used as the parity-check matrix of the code. 

The girth of the ( 7 , p) -biregular bipartite graph Q(V. C) (defined as length of the shortest 
cycle) of the partial geometry PaG( 7 , p, 5 ) for 5 > 1 , is 6 and there are 7777(7 —l )(<5 — l)(p — l )/6 
cycles of length 6 in G(V, C) |f2Tj. For S = 1 the partial geometry is a generalized quadrangle 
(GQ) lf35l and the corresponding graph will have girth of 8 . A simple modification of the 
argument of ETTl to find the number of cycles of length 8 is as follows. Consider a line L and 
points 77, 77 on L . Consider the intersecting bundles of lines, A (77) and A (77). Each set contains 
(7 — 1) lines other than L . Consider a point v on a line L -> through 77 other than L (i.e., in 
A(n 2 )\{L}) and a line L x through 77 in A(77)\{L}. By the definition of a GQ, there is a unique 
line through v that intersects L\. This is true for each point on the lines in A (77), other than 
77 and each line in A( 77 ) other than L. Hence, there are (7 — l) 2 (p — 1) distinct quadrangles 
in the GQ that contain 77,772 as a side and (£) ways of choosing the two points on a line. The 
argument is repeated for each of the m lines in the geometry and note that each quadrangle is 
repeated four times in this count giving the number of quadrangles in the geometry (cycles of 
length 8 in the Tanner graph) as 

m (P~ !)(7- i) 2 ^) / 4 ‘ 

III. A Special Category of Partial Geometries 

In this section, we present a special category of partial geometries which are specified by 
their line-point adjacency matrices. First, we prove the conditions for a matrix of a specific type 
to be the line-point adjacency matrix of a partial geometry. Some structural properties of this 
category of partial geometries are then developed. In Sections [V] and ED we will present two 
methods for constructing this category of partial geometries based on prime fields and cyclic 
subgroups of prime orders of finite fields, respectively. 
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A. Line-Point Adjacency Matrices 

Consider anmxn matrix H with 0- and 1-entries for which two rows (or two columns) have 
at most one place where both have 1 entries. Such a constraint on rows and columns are called 
the row and column constraint (RC-constraint) |[24|. The matrix H is called an RC-constrained 
matrix. An RC-constrained matrix H is (7, p)-biregular if each of its n columns has 7 1-entries 
and each of its m rows has p 1-entries. The parameters 7 and p are called the column and row 
weights of H. respectively. Suppose the rows of a (7, p) -biregular matrix H are labeled from 
0 to m — 1 and the columns from 0 to n — 1. The 7 rows of H that have 1-entries at the j-th 
position with 0 < j < n are said to be attached to the j-th column. Likewise, the p columns 
that have 1-entries in the z-th row are said to be attached to the z-th row. 

The following theorem gives the conditions that an RC-constrained m x n regular matrix H 
with 0- and 1-entries is the line-point adjacency matrix of a partial geometry with n points and 
m lines. 

Theorem 1 . Let H be an m x ?z RC-constrained (7, p)-biregular matrix which is a 7 x p array of 
7x7 circulant permutation matrices (CPMs), where m — y 2 and n = 7 p. Then, H is the line- 
point adjacency matrix of a partial geometry PaG( 7 , p, p — 1) that has n points corresponding 
to the columns of H and m lines corresponding to the rows of H. 

Proof: Since H is a 7 x p array of 7 x 7 CPMs, it is a y 2 x 7 p (7, p)-biregular matrix with 
constant column weight 7 and constant row weight p. It is clear that each point (corresponding 
to a column) is on 7 lines (corresponding to 7 rows) and each line (corresponding to a row) 
passes through p points. The RC-constraint implies that any two points are on at most one line. 

It remains to show that if a point v is not on a line L, then there are exactly p — 1 lines that 
pass through v and p — 1 points on L. We group the n points into p sets of size 7 where the 
points in each set correspond to 7 consecutive columns in H, namely, the columns comprising a 
column-block of CPMs in H. Since every row has p ones, then by adding all the rows attached 
to the column corresponding to the point v, where the sum is over the integers rather than over 
GF(2), we obtain a vector, z, of length n whose components as integers add up to 77. Notice 
that the entry in the column corresponding to the point v in z is 7 while all other entries in 
columns corresponding to the 7 — 1 other points in the same set as v are zeros since every 
row attached to the column corresponding to v has zeros in all columns corresponding to other 
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points in the same set as v as every CPM has a single 1 in each row. Hence, all the (p — 1)7 
components in z other than those corresponding to v and the 7 — 1 points in the same set as v 
add up to (p — 1 ) 7 . Because of the RC-constraint, all these components are at most equal to 1 
implying that all of them are equal to 1. Hence, every column other than those corresponding to 
v or other points in its set is attached to a unique row corresponding to a line passing through v 
while every column corresponding to a point other than v but in the same set as v is not attached 
to any row corresponding to a line passing through v. Since L passes through a single point in 
each set, it passes through p — 1 points in sets other than that of v. Each of these points is on a 
line passing through v. The point on L that is in the same set as v is not adjacent to v. Hence, 
there are exactly p — 1 lines that pass through v and p — 1 points on L. This proves that H is 
the adjacency matrix of a partial geometry PaG( 7 , p, p — 1). ■ 

The proof of the above theorem gives insight into the structure of the partial geometry 
PaG( 7 , p, p — 1) which is discussed further in following sections. 

It is of interest to develop characterizations of partial geometries in the sense that if a 
combinatorial construction is developed, what technique could be used to verify that it is a 
partial geometry? The next discussion considers one such approach for a particular case where 
the adjacency matrix of the configuration is an array of CPMs. 

Let H be the m x n line-point adjacency matrix of a partial geometry PaG( 7 , p, 5 ). Suppose 
m = rn'(i and n = n'£ and assume that H is an m! x n' array of CPMs of order l, i.e., if P is 
a generator of the cyclic matrix group with the first (or the top) row the (-tuple (0,1, 0, • • • , 0). 
The group of such £ x £ matrices is {I, P, ..., P^ 1 }, P' = I, where I is an £ x £ identity matrix. 

The matrix HH T is an m x m matrix of 0’s and l’s whose th element is 1 if the 7-th 
line intersects the 7 -th line and 0 if the lines do not intersect and with diagonal elements p. 
If the i-th column-block of CPMs of H is multiplied by the £ x £ matrix I + P J ’, 0 < i < n' 
and 0 < j < £, the effect is to add to the circulant in a given row of circulants, the power 
of a nonunity permutation (circulant), i.e., it adds a 1 to each row of the circulant in the 7 -th 
column-block of circulants, not in the position of the existing 1. For 0 < i < n' and 0 < j < £, 
form the following n' x n' diagonal array of circulants of size £ x £: 

D^ = diag(I.I,...,I + PV..,I) 



where the matrix I + P' is in the i-th position of the array D' : 'Then, the matrix 

G (i ’ f) = HD (!j) H t , 0 < i < ri, 0 < j < t (3) 

is an m x m matrix with 5 2’s per row, corresponding to the <5 lines that intersect the line of that 
row with the extra point not on the line, diagonal elements of p and all other elements 0 or 1 . 
Note that the matrix contains an extra point in each row, each row corresponding to a 

line and a point not on the line. Conversely, if H is an m x n binary matrix which is an m! x n' 
array of CPMs of order ^ satisfying these conditions, it represents a partial geometry PaG( 7 , p, 5). 
This is clear by construction. Thus, we have the following theorem that characterizes a binary 
array of CPMs to be a line-point adjacency matrix of a partial geometry. 

Theorem 2. Let H be an m x n binary matrix written as an m! x n' array of CPMs of order 
i whose associated matrices G (l%1> of (0) satisfy the stated conditions. Then, H represents a 
partial geometry PaG(y, p, 5). 

The condition of the theorem does not allow zero matrix (ZM) of order i. The argument can 
be modified to accommodate this case although the result is more complicated and not presented 
here. 

In Sections [V] and ED we will use Theorems Q] and [2] to construct two new classes of partial 
geometries whose line-point adjacency matrices are arrays of CPMs. 

B. Structural Properties 

Consider a partial geometry PaG( 7 ,p, p — 1) as in the above Theorem Q] Notice that the 
line-point adjacency matrix H consists of 7 row-blocks of CPMs and p column-blocks of CPMs 
of size 7 x 7 . We first group the 7 p points into p sets, denoted by N 0 , N 1; ..., N p _ l5 each 
consisting of 7 points which correspond to 7 columns in a column-block of CPMs of the line- 
point adjacency matrix H. Since each row of a CPM has only one 1-entry, the p points of a line 
L are distributed in p different sets, one point in each set. For 0 < k < p, a point v k in the set 
Nj that is not on a line L in PaG( 7 , p, p — 1) is not adjacent to the point on L in the same set 
Nj but is adjacent to each of the other p — 1 points on L that are in the p — 1 sets other than 

Nj- 

The y 2 lines corresponding to the 7 2 rows of the line-point adjacency matrix H can be grouped 
into 7 bundles, denoted by M 0 , M 1; ..., M 7 _ l5 each consisting of 7 lines which correspond to 
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the 7 rows in a row-block of CPMs in H. Since any two rows in a row-block of CPMs in H do 
not have any position where they both have 1-entries, the 7 lines in each bundle M, : are parallel 
lines. Hence, each bundle Mj corresponds to a parallel bundle of lines in PaGfy, p, p — 1). 

Let H 0 , Hi,.... H 7 _i denote the 7 row-blocks of H. For 0 < i < 7 and 0 < k < 7 , let 
hj k = (hj fc 0 , ..., tP _i) be the k-\h row in the i-th row-block H, which consists of p 

sections, hj )fe)0 , ..., hj !fcj/9 _i, each consisting of 7 components. Each section is the k -th row 
of a CPM in Hj. If we cyclically shift all the p sections of h, /,. simultaneously one place to the 
right within the sections, we obtain the (k + l)-th row hj fc+ i of Hj which also consists of p 
sections and each section is the (k + l)-th row of a CPM in Hj. The above cyclic-shift within 
each section of h,./, is referred to as section-wise cyclic-shift of the row h,./,. For k — 7 — 1, the 
section-wise cyclic-shift of hj i7 _ x results in the top row hj 0 of Hj. Consequently, all the rows 
of the i-th row-block Hj can be obtained by section-wise cyclically shifting the top row h,,n of 
Hj 7 — 1 times. 

Let Hg be an 7 x 7 p matrix over GF(2) which consists of the top rows h 0 ,o, h 10 ,..., h 7 _ 10 
of the 7 row-blocks H 0 , Hi,..., H 7 _i of the line-point adjacency matrix H of PaG(y, p, p — 1). 
Then, it follows from the section-wise cyclic structure of H, the entire array H can be obtained 
by section-wise cyclically shifting Hq 7 — 1 times. This structural property is referred to as QC 
(or section-wise cyclic) structure. The QC-structure of H implies that the incidence vectors of 
the y 2 lines in the partial geometry PaG( 7 , p, p — 1) which correspond to the y 2 rows of H can 
be generated by the incidence vectors of the 7 lines which correspond to the 7 rows in Hq. That 
is to say that the lines corresponding to the rows of Hq and the points on these lines specify 
the partial geometry PaG(y, p, p — 1). 

Due to the QC-structure of H, we call the partial geometry PaG(y,p, p — 1) a quasi-cyclic 
partial geometry, denoted by QC-PaG(y,p, p — 1). 

In the following, we give an example to demonstrate that a QC-PaG does exist. 


Example 1. Consider the 3x3 array H of CPMs of size 3x3 given in the equation below: 
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H = 


1 0 0 
0 1 0 
0 0 1 
1 0 0 
0 1 0 
0 0 1 
1 0 0 
0 1 0 
0 0 1 


1 0 0 
0 1 0 
0 0 1 
0 1 0 
0 0 1 
1 0 0 
0 0 1 
1 0 0 
0 1 0 


1 0 0 
0 1 0 
0 0 1 
0 0 1 
1 0 0 
0 1 0 
0 1 0 
0 0 1 
1 0 0 


(4) 


It is a 9x9 matrix over GF(2). By checking, we find that it satisfies the RC-constraint. Hence, 
it follows from Theorem [7] that H is a line-point adjacency matrix of a PaG(3, 3, 2) with 9 points, 
9 lines and connection number 2. Each line consists of 3 points and each point is on 3 lines. 
The lines in PaG( 3, 3, 2) can be partitioned into 3 parallel bundles, each consisting of 3 parallel 
lines. Suppose we use Vo,vi,v 2 ,v 3 ,V 4 ,V 5 ,ve,v 7 ,vs to denote the 9 points in PaG( 3,3,2) which 
correspond to the 9 columns of H. Then, the 9 lines in PaG( 3,3,2) are: 


L 0 = {v 0 , v 3 , v 6 }, 
L\ ={vi,V 4 ,V 7 }, 

L 2 = {v 2 ,v 5 ,v s }, 
L 3 ={vo,V4,Vs}, 
U ={vi,v 5 ,v 6 }, 
L 5 = {v 2 ,v 3 ,v 7 }, 
Lq ={v 0 ,v 5 ,v 7 }, 
L 7 = {Vi,V 3 ,U 8 }, 
Lg ={v 2 ,v 4 ,v 6 }. 


Group the points in PaG(3, 3, 2) into 3 disjoint sets, N 0 ,Ni and N 2 , each consisting of 3 
points. For 0 < j < 3, the 3 points in N ? correspond to the 3 columns in the j-th column- 
block of H. Then, N 0 = {w 0 , v\, v 2 }, N x = {v 3 , v 4 , v 5 } and N 2 = {vq,v 7 ,v$}. Consider the line 
L.> = {t; 0 , v 4 , t' 8 }. The point v 5 is not on L 3 . By checking, we find that this point is adjacent to 
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the points vq and v%. Notice that the point v§ is in the set Ni but its two adjacent points vq and 
v 8 are in two separate sets N 0 and N 2 . The graphical representation of PaG( 3,3,2) is shown 
in Fig. [7] and the Tanner graph associated to the line-point adjacency matrix H is shown in Fig. 

1 2} The girth of the graph is 6. 



Fig. 1. A graphical representation of the partial geometry PaG(3, 3, 2) specified by its line-point adjacency matrix given by 0 
in Example [I] 


Vo Vi V 2 V 3 V A V 5 V 6 V 7 V 8 



Fig. 2. The Tanner graph associated with the line-point adjacency matrix of PaG(3, 3, 2) given by 0 in Example Q] 


IV. Protograph Representation of a QC-PaG 

As shown in Section HU a partial geometry with n points and m lines can be represented by 
a bipartite graph with n VNs which represent the n points of the geometry and rn CNs which 
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represent the m lines of the geometry. A VN v 3 is connected to a CN c, by an edge (i,j) if and 
only if the point represented by the VN v 3 is on the line represented by the CN c r . For a partial 
geometry with a large set of points and a large set of lines, the bipartite graph representation of 
the geometry would be very large and complex. However, the QC-PaG( 7 , p, p — 1) characterized 
by Theorem 1 can be effectively represented by a much smaller bipartite graph based on its 
QC-structure. 

Recall that the line-point adjacency matrix H of the QC-PaG( 7 ,p, p — 1) is a binary 7 x p 
array of CPMs of order 7 of the following form: 


H 


Ho 

Hi 


H 7 _i 


> 

0 

0 

O 

< 

> 

O 

1 

Ai,o 

A 1(1 • • 

Ai )P _i 

> 

1 

O 

A 7 _i ; i • • 

> 

1 

1 


(5) 


The array H = [-^i,j]o<i<'y,o<j< P consists of p column-blocks of CPMs, denoted by N 0 , ..., 

N p _i, and 7 row-blocks of CPMs, denoted by H 0 . H 1; .... H 7 _!. For 0 < i < 7 and 0 < j < p, 
each column-block Nj consists of 7 consecutive columns of H and each row-block FT, consists 
of 7 consecutive rows of H. The 7 columns of the j-th column-block Nj correspond to 7 points 
of the QC-PaG( 7 , p, p — 1) and the 7 rows of the z-th row-block H correspond to 7 lines of 
the QC-PaG( 7 , p, p — 1). 

For 0 < j < p, let ( l> 3 be the set of 7 points which correspond to the 7 columns of the j-th 
column-block Nj of H. For 0 < i < 7 , let 0, : be the set of lines which correspond to the 7 rows 
of the z-th row-block H, of H. In forming the Tanner graph <j PaG of the QC-PaGfy , p, p — 1), 
the 7 points in <I> ? are represented by 7 VNs and the 7 lines in 1 1, are represented by 7 CNs. 
Hereafter, we use points and VNs interchangeably, lines and CNs interchangeably. The 7 VNs 
in <f>j are called type-j VNs and the 7 CNs in Q, are called type-/’ CNs. From the QC-structure 
of H, we see that a type-j VN can only be connected to a type-/ CN and vise versa. In forming 
the Tanner graph <j PaG of the QC-PaG( 7 , p, p — 1), the 7 type-j VNs in <3> ? are connected to the 
7 type-i CNs in 0, : based on the 7 1-entires in the CPM A* j and vise versa. If we label the 
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columns and rows of a CPM A ? J in H from 0 to 7 — 1, then A it j is uniquely specified by the 
location of the single 1-entry of its top row, called the generator. If the single 1-entry of the top 
row of A ij locates at the position k tJ , 0 < k, 3 < 7, then we use (k l3 ) to specify the CPM A h3 . 

With all the terms defined above, we now construct a bipartite graph, denoted by f/ protoPaG , 
with p VNs and 7 CNs as shown in Fig. 0 The p VNs of (/ proto , p a o represent the p clus¬ 
ters $ 0 , $ 1 , • • •, $p-i of VNs in Cp a (, and the 7 CNs in ty pro to,PaG represent the 7 clusters 
O 0 , fli,..., fl 7 _i of CNs in Ci'aG• In Cpmto.PaG, the VN is connected to the CN Q, by an 
edge labeled by (kif) which is the location of the single 1 -entry of the generator (or top row) 
of the CPM A ij in the line-point adjacency matrix H of the QC-PaG( 7 , p, p — 1). The bipartite 
graph f?p r oto,PaG contains all the structural information of the QC-PaGfy , p. p — 1). The size of 
the bipartite graph Oproto.PaG is smaller than the size of the bipartite graph 0 p aG by a factor 7 . 
This bipartite graph Opmto.PaG is called the protograph of the QC-PaG( 7 , p, p - 1). Basically, the 
protograph f/ pro to,PaG of the QC-PaG( 7 , p, p — 1) consists of p super-VNs, $ 0 , $ 1 , • • •, $ P -i, and 
7 super-CNs, O 0 , Q 1 ,..., 0 7 _p. The Tanner graph C/ PaG of the QC-PaG( 7 , p, p — 1) is simply an 
expansion of the protograph (? pro to.PaG of the QC-PaG( 7 ,p, p — 1). 


$0 $1 ® P -1 



Fig. 3. The protograph of the QC-PaG( 7 , p, p — 1). 


Example 2. Consider the QC-l J aG(3. 3, 2) given in Example 0 We group the 9 points of the 
geometry into three disjoint sets, <T 0 = {v 0 ,Vi,v 2 }, = {73, rq, v 5 }, $ 2 = {^6, ^7, r> 8 }, and 

the 9 lines of the geometry into three disjoint sets, flo = {Lq, L 2 }, Oi = {L 3 , L 4 , L-,}, 
= {L$, L 7 , L 8 }. Using <h 0 > ^ 1 , < I ) 2 cis the super-VNs and f2 0 > ^2 cis the super-CNs, 

we form the protograph Q pro to,PaG of the QC-PaG( 3, 3, 2) as shown in Fig. 0 Each edge is 
labeled by the position of the single 1-entry of the generator of each CPM which specifies the 
connections between the constituent VNs of a super-VN ( l>j in Q pro to,PaG a nd the constituent CNs 
of a super-CN Q, in Gpmto.Pac■ 
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$0 $ 1 $2 



Fig. 4. The protograph of the QC-PaG(3,3, 2) given in Example [2] 

LDPC codes can be constructed based on relatively small protographs Q], lfl3ll . Ifl4l . fl47ll . 
These LDPC codes are called protograph LDPC codes. In constructing such a code, we take 7 
copies of a chosen protograph. Then, we permute the edges of the individual copies and connect 
the copies into a large Tanner graph. The null space of the adjacency matrix of the Tanner 
graph gives a protograph LDPC code. Constructions of protograph LDPC codes in most of the 
literature are pseudo-random and require computer aid to perform the connections of the copies 
of the chosen protograph. 

From the graph point of view, the QC-LDPC code constructed based on the QC-PaG( 7 , p. p - 
1) may be regarded as a protograph QC-LDPC code. Based on the description of the QC- 
PaG( 7 ,p, p— 1), we first construct its protograph C pr oto,PaG- Then, we take 7 copies of f? pro to,PaG, 
permute the edges of the copies and connect the copies to form the Tanner graph Civc, of the 
QC-PaG( 7 ,p, p — 1). The null space of the adjacency matrix of the Tanner graph Civr, of the 
QC-PaG( 7 , p, p — 1) gives a protograph QC-LDPC code. In Sections lYl and IVTl we will present 
two classes of QC-PaGs. The protograph of a QC-PaG from each of these two classes can be 
constructed directly from the description of its base matrix. 

V. A Class of QC-PaGs Constructed Based on Prime Fields 

In this section, we present a class of QC-PaGs based on prime fields whose line-point adjacency 
matrices are arrays of CPMs. 

A. Construction 

Let p be a prime and GF(p) be a prime field which consists of the following p elements: 
0,1, 2 ,...,p — 1. We represent each element i in GF(p) by a p x p CPM, denoted by Q(*), with 
columns and rows labeled from 0 to p — 1 , whose generator has its single 1 -component at the 





15 


location i. This representation is one-to-one. For i — 0, Q(0) is a p x p identity matrix. Q (i) is 
referred to as the CPM-dispersion of the element i in GF(p). 

Form the following p x p matrix B p over GFf/i) with columns and rows labeled from 0 to 

p — 1: 



0-0 

0-1 

0 • (p 

-1) 


1-0 

i-i 

1 • (p 

-1) 

B p = 

2-0 

2-1 

2 • (p 

-1) 


1 

o 

('P - 1) ■ 1 • 

■ (P ~ 1) • 

(P~ 


i.e., B p = \bij]o<i,j<p- i with b t;i = ij, where the multiplication of two elements in GF(p) is 
carried out under modulo p. 

The matrix B p has the following structural properties: 

1) All the entries in the 0-th row and column are zeros; 

2) All the entries in any row or column, other than the 0-th row and column, are different 
and contain all the elements of GF(p); 

3) Two different rows or columns have the 0 element of GF(p) in common at the 0-th row 
or column and differ in all the other pm 1 positions; 

4) The ?’-th column of B p is identical to the transpose of i-th row of B p for 0 < i < p. 

The last property implies that the transpose Bj of B p is identical to B p , i.e., B P T = B p . 

If each entry in B p is dispersed into its corresponding pxp CPM, we obtain a p x p array H p 
of CPMs of size p x p. It is a p 2 x p 2 matrix over GF(2). Each row has weight p and each column 
has weight p. Based on the structural properties of B p , we can readily see that H p satisfies the 
RC-constraint. Hence, it follows from Theorem Q] that H p is the line-point adjacency matrix of 
a QC-PaG P (p,p,p — 1) with columns corresponding to the points and rows corresponding to 
the lines in the QC-PaG p (p, p, p — 1 ). Equally, Theorem [ 2 ] could have been applied. The partial 
geometry QC-PaG P {p,p,p — 1) consists of n = p 2 points and m = p 2 lines. Each line in the 
QC-PaG p (p,p, p — 1) consists of p points and each point is on p lines. A point v that is not on 
a line L is connected to p — 1 points on L by lines, i.e., the connection number of a point in 
the QC-PaG P {p,p,p — 1) is p — 1. The array H p is called the CPM-dispersion of B p and B p is 
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called the base matrix for the construction of the QC-PaG p (p, p, p — 1). The subscript “p” of the 
QC-PaG P (p,p,p — 1) stands for “prime”. 

As an array of CPMs, H p consists of p row-blocks of CPMs and p column-block of CPMs. 
Two different rows in a row-block of H p have no position where they both have 1-entries. 
Hence, the p rows in each row-block correspond to p parallel lines in the QC-PaG P (p,p, p — 1). 
Therefore, the QC-PaG P (p, p,p — 1) consists of p parallel bundles, each consisting of p parallel 
lines. Let M 0 , Mi,..., M p _i denote the p parallel bundles of lines in the QC-PaG P (p,p,p — 1). 
For 0 < i < p, the p lines in the /-th parallel bundle M, correspond to the p rows in the 
i-th row-block of H p . The p lines in each parallel bundle contain all the p 2 points in the QC- 
PaG P (p,p,p — 1). This partial geometry is a net. 

The p rows in H p that have 1-entries at the j-th position are in p different row-blocks. 
These rows correspond to p lines in an intersecting bundle of lines that intersect at the point Vj 
corresponding to the j-th column of H p . The p lines in an intersecting bundle are in p different 
parallel bundles of the QC-PaG P (p,p,p — 1). 

The p 2 points in the QC-PaG P (p,p,p— 1) can be divided into p sets, N 0 , Nj,..., N p _i, each 
consisting of p points. For 0 < j < p, the p points in the j-th set N ? correspond to the p columns 
of the j-th column-block of H p . The p points of a line L are distributed in p different sets, one 
and only one point in each set. A point v in the set N y that is not on a line L is adjacent to 
one and only one point on L in each set other than Np 

It follows from the results on the structural properties given in [ITHl that the Tanner graph, 
denoted by C? P)PaG , of the QC-PaG P (p,p,p — 1) has girth 6 and contains p 3 (p — l) 2 (p — 2)/6 
cycles of length 6. 

The protograph, denoted by f? PiPro to,PaG of the QC-PaG p (p,p,p — 1) can be constructed directly 
from the base matrix B p given by ©. It consists of p VNs and p CNs, labeled by <f> 0 , Tq, ..., $ p _! 
and O 0 , , Q p _i, respectively. The p VNs and p CNs of the protograph <j PiPro to,PaG correspond 

to the p columns and p rows of the base matrix B p of the QC-PaG p (p,p, p — 1). Every VN is 
connected to every CN. The edge (j, i) connecting the VN ( l> 3 to the CN Q, is labeled by the 
integer ij modulo p. 

If we remove the first column and the first row from the base matrix B p , we obtain a Latin 
square of order p — 1. The use of Latin squares for constructing LDPC codes was considered in 

m, G33. 
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Example 3. Let p = 3 and let GF( 3) be the field to construct a QC-PaG p { 3, 3, 2). Using (0), we 
find the base matrix for constructing the line-point adjacency matrix H p of the QC-PaG v {?>. 3, 2) 


is 


The 3 x 3 CPM-dispersion H p of B p is 


H r 





0 

0 

0 




B 

p = 


0 

1 

2 







0 

2 

1 




B p 

is 








1 

0 

0 

1 

0 

0 

1 

0 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

0 

1 

0 

0 

1 

0 

0 

1 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

1 

0 

0 

0 

1 

1 

0 

0 

0 

0 

1 

1 

0 

0 

0 

1 

0 

1 

0 

0 

0 

0 

1 

0 

1 

0 

0 

1 

0 

1 

0 

0 

0 

0 

1 

0 

0 

1 

0 

1 

0 

1 

0 

0 


(7) 


( 8 ) 


which is exactly the line-point adjacency matrix of the PaG( 3,3,2) give?? £>y (0) in Example\2\ 
The Tanner graph and the protograph of this partial geometry are shown in Fig. | 2] and Fig. [7] 
respectively. 


B. Subgeometries 

Let r be a positive integer less than p — 1. Suppose we delete r column-blocks from H ; ,. We 
obtain a p x (p — r) subarray, denoted by H p (p,p — r) , of H p . It follows from the structural 
property developed above that H p (p, p—r) is the line-point adjacency matrix of a QC-PaG p (p, p— 
r, p — t — 1) which has p(p — r) points and p 2 lines, each point on p lines and each line consisting 
of p — t points. The connection number of a point in the QC-PaG p (p, p — r, p — r — 1) is p — r — 1. 
The QC-PaG p (p, p — r, p — r — 1) is a subgeometry of the QC-PaG p (p, p,p — l). Its Tanner graph 
(or protograph) is a subgraph of the Tanner graph (or protograph) of the QC-PaG P {p,p,p — 1). 

Therefore, for a given prime field, we can construct a family of QC-PaGs. 
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C. QC-LDPC Codes on the QC-PaG p (p,p,p — 1) 

For 1 < k,r < p, let B p {k,r) be a k x r submatrix of the matrix B p of ©. CPM-dispersing 
the entries of B p {k,r), we obtain a k x r array H p (fc,r) of CPMs of size p x p which is a 
subarray of H p . The array H p (k,r) is a kp x rp matrix with column and row weights, k and 
r, respectively. The null space of H p (k,r) gives a (k, r) -biregular QC-LDPC code, denoted 
by C P!qc (k,r), of length rp with rate at least (r — k)/r. Therefore, for a given prime field, a 
family of QC-LDPC codes of various lengths and rates can be constructed. The Tanner graph 
of C P}qc (k, r ) has girth of at least 6. 

Express B p (k,r) as B p (k,r) = [bij]o<i<k,o<j<r- Suppose an entry in B p (k,r) is replaced 
by the zero element of GF(p). In the CPM-dispersion H p (k,r) of B p (k,r), this replacement 
results in replacing a p x p CPM in H p (A\ r) by a p x p zero matrix (ZM) which is referred to 
as masking |fTTfl. Ifl2l . |j42fl. [f53l Let A be a nonnegative integer less than the number of total 
nonzero entries in B p (A\ r). The replacement of A nonzero entries in B p (k, r ) by A zeros amounts 
to replacing A CPMs by A ZMs at the locations in H p (k,r) corresponding to the locations of 
the A entries in B p (fc,r) which are replaced by zeros. Masking A CPMs in H p (fc,r) amounts to 
removing A p edges from the Tanner graph Q p (k,r) associated with H p (k,r). Removing these 
edges in Q p (k,r) may break many short cycles in Q p {k. : r). As a result, the resultant Tanner 
graph f? Pi mask(fc, r) may have a much smaller number of short cycles, or a larger girth, or both. 
The subscript “mask” stands for “masking”. In choosing the entries in B p (fc,r) to be masked, 
we have to avoid disconnecting the Tanner graph of H p (A:, r). 

The operation of masking B p (A;,r) = |Ajj]o<i<a-.o</. ;r can be modeled mathematically. Let 

Z(k. r) = [£ij]o<i<fc,o<j<r be a k x r matrix with the zero element and unit element of GF(p) 
as entries. Define the following product of Z(k,r) and B p (A;,r) : B Pjinask (A :,r) = Z(k,r) © 

B p (k,r) = [©,/A ; j-]o<i</,-.o</ :r (the Hadamard product) where Zijbij = b i:j if Zij = 1 and 

z i,jbi,j = 0 if Zij = 0. In this matrix product operation, entries in B ; ,(A\ r) at the locations 
corresponding to the locations of zero-entries in Z(A:, r) are replaced (or masked) by 0’s. The 
CPM-dispersion of B Pimask (A;, r) gives a k x r masked array H P)mask (A:, r) of CPMs and ZMs 
of size p x p. We call Z(k. r) and B p mask (/c,r) the masking matrix and the masked base 
matrix, respectively. The null space of H pmask (A:, r) also gives a QC-FDPC code, denoted by 
G Pi(?Cimask (A;, v). 

Example 4. Consider the 127 x 127 base matrix B p over GF{ 127) given in It contains the 
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following 4x8 submatrix B p (4, 8) (rows and columns chosen at random): 


2 

83 

33 

46 

36 

94 

42 

86 

109 

15 

84 

94 

57 

43 

3 

115 

112 

76 

70 

36 

111 

57 

66 

117 

31 

80 

67 

78 

50 

60 

16 

63 


We design the following masking matrix: 

0 1 1 1 l" 

11111 
110 10 
1 0 1 0 1_ 

Masking B p (4, 8) with Z(4, 8) gives the matrix B p ma ^(4, 8). Replacing each nonzero entry 
in B p ma ^(4, 8) by its corresponding 127 x 127 CPM and each O-entry by a 127 x 127 ZM, we 
obtain a masked 4x8 array H p>maifc (4, 8). It is a 508 x 1016 matrix with column and row weights 
3 and 6, respectively. The null space o/H p masfc (4, 8) gives a (3, 6) -biregular binary (1016,508) 
QC-LDPC code with rate 0.5. The Tanner graph of this code has girth 8 and contains 889 
cycles of length 8. The bit and block error performances of the code decoded with 5, 10 and 
50 iterations of the min-sum algorithm (MSA) l \TU\I are shown in Fig. |5| We see that the code 
has a very low error-floor. Included in the same figure are the bit and block performances of 
an LDPC code constructed using the progressive edge growth (PEG) algorithm fdV\l . The PEG 
code is only decoded with 50 iterations of the MSA. We see that the code constructed based on 
GF( 127) outperforms the PEG code below the bit error rate (BER) of 10 “ 7 . 

As noted, a special case of the partial geometries is when 5 = 1, the generalized quadrangles 
(GQs). In this case there are no triangles in the geometry and hence the associated graph has 
girth 8. There have been numerous studies on the use of GQs in coding theory (although none 
considered their trapping sets to our knowledge). The works include J2]-[[4|, lf30l . 01 . Il36ll . 
m, ffl, EEL El. Some comments on the codes from GQs derived from these papers are 
given. The paper of Liu et al 01 contains a very useful summary of the minimum distance 
bounds and rates of the various codes resulting from using GQs for code constructions. Among 
other things, that work established the following two important results: 

1) A code whose parity-check matrix is the adjacency matrix of a generalized rl-gon, for n 


Z(4, 8) = 


1 0 1 
0 1 0 
111 
111 
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Fig. 5. Performance curves for the (1016, 508) QC-LDPC code and the PEG code given in Example [4] 


even, has only codewords of even weight. 

2) A code whose parity check-matrix is the adjacency matrix of a GQ is quasi-cyclic. 

This last point is important for code implementation. Its proof considers the transitivity of the 
automorphism group of the quadrangle which produces different sized circulant blocks for the 
different constructions. The reader is referred to m for details. 

VI. A Class of QC-PaGs Constructed Based On Cyclic Subgroups of Finite 

Fields 

In this section, we present another class of QC-PaGs which are constructed based on the 
cyclic subgroups of finite fields. 

A. Construction 

Let GF(g) be finite field with q elements where q is a power of a prime. Let t be a prime 
factor of q — 1 and q — 1 = ct. If a is a primitive element of GF(g), then (3 = a c is an element of 
order t in GF(g). The set S = {l,/3,/3 2 ,... ,(3^^} forms a cyclic subgroup of GF(g) of order 

t. 
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Form the following t x t matrix: 


B c = 


1 1 

1 p 

1 p 2 


1 

P 2 

C 8 2 ) 2 


1 


pt-i 

m - 1 


(9) 


1 / 3 t_1 (/ 3 4 - 1 ) 2 ••• (/ 3 t - 1 ) f - 1 


Express the above matrix in the form B c = [P^]o<i,j<t whose (z.j)-th element is p i P The 
subscript “c” stands for “cyclic group”. All the entries of B c are elements of the cyclic subgroup 
S. Let P denote the CPM of size txt whose first (or the top) row is the t-tuple (0,1, 0, ■ • • , 0) 
with the single 1-component at the position-1 (the positions are labeled from 0 to t — 1). For 
0 < i < t, the Ath power of P. denoted by P f is a CPM of size txt whose top row has its 
single 1-component at the position-^. Then, {I, P,..., P (t_1 )} is a cyclic matrix group of order 
t with P as a generator and P 4 = P° = I. 

For 0 < i < t, we represent the element p l in the cyclic group S of GF(g) by the CPM P . 
This representation is one-to-one and P p is the CPM-dispersion of d ( '. 

Let H c = B C (P) denote the CPM-dispersion of the matrix B c , i.e., for 0 < i,j < t, the 
entry /T- 7 of B c at the location (i,j) is dispersed into the CPM P ; where i = (ij) modulo t. 
Then, H c = B C (P) is a txt array of CPMs of size txt. This array H c = B C (P) is a t 2 x t 2 
binary matrix which is the line-point adjacency matrix of the QC-PaG c (t, t, t — 1). The result 
can be established using either Theorem Q] or Theorem [2l but Theorem [2] is used for illustrative 
purposes. Notice that P 7 = P _1 . Consider the product of B C (P) and its transpose [B C (P)] T . 
Express the product B C (P) [B C (P)] T as a txt array of t x t matrices. Then, we have 


A = [A r , fl ] 0 < r , a < t = B c (P)[B c (P)] t . 


For 0 < r,s < t, the constituent matrix A r s of A is 


A 

-^rjS 


= eI; 1 oP"(pT s 



if r s 
if r = s 


where J is the t x t all ones matrix and I the txt identity matrix. Suppose the ?-th column-block 




22 


of CPMs of B(P) is multiplied by I + P J to give A^ lJ) = (A ' s ) where 

a; >s = A r , s + p i(r - s) +L 

Thus, each CPM of each row-block of B C (P) has a CPM added to it. For a non-diagonal 
circulant (r ^ s ), this means each row has (t — 1) 2’s with all remaining elements in the row 
being 0 or 1. It follows from Theorem [2] that the array H c = B C (P), as a t 2 x t 2 matrix, is 
the line-point adjacency matrix of a partial geometry PaG c (t, t, t — 1). The PaG c (f, t, t — 1) is a 
QC-PaG consisting of n = t 2 points and m = t 2 lines with connection number 5 — t — 1. This 
partial geometry PaG c (t,t,t — 1) is also a net. 

If we want to use Theorem Q]to prove H c = B C (P) is the line-point adjacency matrix of the 
partial geometry PaG c (f, t, t — 1), we need to show that any 2x2 submatrix of B c is nonsigular. 
This can be done easily. It follows from lf55l Corollary 1] that H c as a t 2 x t 2 matrix satisfies 
the RC-constraint. Consequently, it follows from Theorem Q] that H c = B C (P) is the line-point 
adjacency matrix of the partial geometry PaG c (t,t,t — 1). 

The above construction gives another class of QC-PaGs whose line-point adjacency matrices 
are arrays of CPMs. For 0 < r < t, if we delete r column-blocks from H c = B C (P), we obtain 
a t x (t — t) array H c (f, t — r) of CPMs of size t x t which is a t 2 x t(t — r) matrix. H c (f, t — r) 
is the line-point adjacency matrix of a partial geometry PaG c (t,t — r, t — r — 1) which is a 
subgeometry of the PaG c (t,t,t — 1). 

Notice that all the entries in any row or column, other than the 0-th row and column, are 
different, since the order of (3 is a prime. By removing the 0-th row and column from B c , we 
obtain a Latin square of order t — 1. 

B. QC-LDPC Codes on the QC-PaG c (t,t,t — 1) 

For 1 < k,r < t, let B c (fc,r) be a k x r submatrix of the matrix B c of (9). CPM-dispersing 
the entries of B c (fc, r), we obtain a k x r array H c (k,r) of CPMs of size t x t which is a 
subarray of H c . The array H c (fc, r) is a kt x rt matrix with column and row weights, k and 
r, respectively. The null space of H < : (k,r) gives a (fc, r) -biregular QC-LDPC code, denoted by 
C c ,q C (k, r), of length rt with rate at least (r — k)/r. The Tanner graph of C Ctqc (k, r ) has girth of 
at least 6. 

In the following, we give an example to illustrate the construction of a QC-LDPC code with 
QC-structure based on the QC-PaG — 1). 
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Example 5. Consider the field GF(2 7 ). Since 2 7 — 1 = 127 is a prime, it cannot be factored. So, 
we choose t = 127. Let a be a primitive element ofGF(2 7 ). Then, the set S = {1, a, a 2 ,..., a 126 } 
forms the only cyclic subgroup of GF{ 127). The order of S is 127. With the choice of t = 127, 
we construct a 127 x 127 base matrix B c = Q iJ ]o<i,?<i 27 in the form of ©. All the 127 nonzero 
elements in GF(2 8 ) appear in each row and each column of B c . Dispersing each entry in B c 
into a 127 x 127 CPM, we obtain a 127 x 127 array H c = B C (P) of CPMs of size 127 x 127. 
The array H c is a 16129 x 16129 matrix with both column and row weights 127. This matrix 
is the line-point adjacency matrix of a QC-PaG r (127. 127,126) with 16,129 points and 16,129 
lines with connection number 126. 

For 1 < k, r < t, let J2 r (k, r) be a k x r submatrix of the base matrix B c . Then, the null 
space of the CPM-dispersion ofB c (k,r) gives a QC-LDPC code of length 127r whose Tanner 
graph has girth of at least 6. 

Label the rows of the base matrix B c from 0 to 126. Suppose we take row-1 to row-6 from 
B c to form the following 6 x 127 submatrix B c (6,127) of B c ; 


B c (6,127) = 


la a 2 
1 a 2 (a 2 ) 2 
1 a 3 (a 3 ) 2 
1 a 4 (a 4 ) 2 
1 a 5 (a 5 ) 2 
1 a 6 (a 6 ) 2 


a 126 

a 2 ] 

j 126 

a 3 ; 

jl26 

a 4 ' 

j 126 

a 5 ; 

j 126 

a 6 ' 

j 126 


( 10 ) 


Notice that matrix B c (6,127) given by ( 1/(71 ) is actucdly the parity-check matrix of the (127,121, 7) 
Reed-Solomon code of symbol length 127 over GF(2‘) with minimum distance 7. Dispersing each 
entry of B c (6,127) into a 127 x 127 CPM, we obtain a 6 x 127 array H c (6,127) of CPMs of size 
127 x 127 which is a subarray of H c . It is a 762 x 16129 matrix with column and row weights 
6 and 127, respectively. The rank of this matrix is 757. Then, the null space of H c (6,127) gives 
a (6,127 )-biregular (16129,15372) QC-LDPC code C cac { 6,127) with rate 0.953. 

The bit and block error performances of this code decoded with the MSA are shown in Fig. 
[6] (computed with an FPGA decoder). We see that the code achieves a BER of 10“ 15 and a 
block error rate (BLER) of almost 10" 12 without a visible error-floor. It has a beautiful waterfall 
performance. Fig. [6] shows the error performances of the code decoded with 5, 10 and 50 
iterations of the MSA. We see that the decoding of the code converges very fast. At the BER 
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of nr the performance gap between 5 and 50 iterations of the MSA is about 0.5 dB. The 
performance curves of the code decoded with 10 and 50 iterations of the MSA almost overlap 
all the way down to the BER of 10 12 . 



Fig. 6. Performance curves for the (16129,15372) QC-LDPC code given in Example [5] 


VII. Trapping Sets and Finite Gepmetries 

Let Q be the Tanner graph of a binary LDPC code where the code is defined as the null space 
of an m x n binary matrix. It is assumed that the n vertices on the left of Q are VNs of the code 
and the m vertices on the right are CNs. The code is to be decoded via a BP decoder, such as 
the sum-product algorithm (SPA) or the MSA. The behavior of these algorithms depends on the 
cycle structure of the code as the iterations of the algorithm attempt to converge to a codeword. 
A trapping set for the code is defined as follows |[25l . If37l : 

Definition 1. 1) A (k, t) trapping set for the code defined by the Tanner graph Q is a subset 

A, | A | — k, of VNs such that the subgraph of Q induced by the set A, denoted by Q[ A], 
has exactly r odd degree CNs (and an arbitrary number of even degree CNs). 

2) The trapping set is said to be elementary if all the associated CNs in Q[ A] have degree 1 
or 2. 

3) The trapping set is called small if n < y/n (n is the code length) and t/k < 4. 








25 


A good description of the motivation for these definitions is given in li25l . [f37ll and the reader 
is referred there for the discussion. In essence, if there is a small trapping set, there is a relatively 
higher probability that the BP algorithm will fail to converge as it iterates. Beyond trapping sets, 
the notion of an absorbing set was introduced in lfl5ll to aid in the convergence analysis but as 
noted, this notion will not be considered in this work. 

Trapping sets have been discussed for a variety of LDPC codes obtained from combinatorial 
structures (in a variety of ways) (see the works llTDl - lfl2l . Il23ll . Il25l - lf27l . Il29l0 . The question 
of interest is to discuss them in light of the structure of the finite geometric structure of the code 
to determine if sharper bounds can be found over those for other codes. 

A useful general result is given in Theorem 2 of ifTTl . llT2ll . It shows that if Q is the Tanner 
graph of an LDPC code with girth of at least 6, with VNs of degree 7 , and Q contains a (ft, r) 
trapping set and k < 7 , then 

r > (7 + 1 — ft)ft. 

A partial geometry code is the dual space to the row space of the m x n adjacency matrix of 
the partial geometry PaG(y, p, 5). 

The trapping sets of such codes were examined in IfTTl . Ifl2l and a brief discussion of those 
results was given. It is noted that the geometric structure of the partial geometry allows a better 
analysis of the trapping sets over matrices derived from less structured objects. Let A be a 
subset of VNs of size ft and Q[ A] the subgraph of the Tanner graph of a code generated by a 
partial geometry. The properties of the trapping set depends on the structure of this subgraph 
induced by A. Since the CNs generated by A correspond to lines through the points of A, it 
will generate a (ft, r) trapping set if there are precisely r lines in Q[ A] passing through an odd 
number of points of A. Let m, be the number of lines of Q\ A] passing through i points of A 
(which is the number of CNs of degree i). Then: 

Theorem 3. ( /[72| Theorem 3]) Let Gpog be the Tanner graph of a PaG( 7 , p, 5). If IS. is a (ft, t) 
trapping set and k < 7 then: 

r A (7 + 1 — ft) ft + (i — 1 ) 2 m, + i{i — 2)?«* 

i odd i even 

and equality holds if 5 = p and the sums go to 2 [(ft + 1 )/ 2 J — 1 and 2 [ft/ 2 j, respectively. 


Note that the first term on the right hand side of this expression is the general bound noted 
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earlier. Also the number of edges in the graph generated by a given set of n VNs is, by definition 
of the partial geometry, K'y. Hence 

irii t = K'f and rn l = r. 

all i odd i 

In the particular case of a net the bound can be improved somewhat. Let L 0 , Li,..., L p _i be 
a set of parallel lines in PaG( 7 , p, 7 — 1) and for a set A of VNs let k* =| A fl L t |. Then: 

Theorem 4. If the set A of VNs of a net is a (k, t ) trapping set, then 

r > (7 — l)ft — n 2 

+ J2i =1 I : 1 odd } I . 

It is noted that the bound agrees with the previous one whenever Ki < 2 for all i and improves 
on it in other cases. A lower bound on the average size of r is also given in Ifl2l Corollary 1], 

It was observed ll29ll that the structure of the partial geometric code allows comments on 
trapping sets for certain types of sets of VNs and a few of these comments are noted here. 

A partial ovoid of a partial geometry PaG( 7 , p, 5) is a set S of points (VNs) such that every 
line in the geometry is incident with at most one point of S. Such a set can have at most 
1 + (7 — l)(p — 1) points. Thus CNs (lines) corresponding to such a set of size k has ^7 lines 
that intersect the k VNs and each such CN has degree 1. Hence rri\ = nn and all other m[s are 
0. The above bounds are difficult to work with in this case but in the case that t > 3 the ratio 
n/r = 1 + t > A and such a set of VNs cannot be a small trapping set. 

Consider a set of k VNs that are colinear, i.e., one line contains all k points. The subgraph 
generated by such a configuration has r = (7 — 1)k + 1 CNs, m K — 1, mi — 7 — 1) and 

t/k = (7 — 1) + 1 / k. For 7 > 4 this does not correspond to a small trapping set. 

From these arguments it is clear that the smallest trapping sets will arise from CN sets that 
contain large numbers of lines between the r points. An extreme case of this is a clique. A 
clique of a partial geometry is a set of VNs which are mutually colinear. A clique of size k 
in the partial geometry will generate a subgraph with (7 — k)k = mi CNs of degree 1 and 
( 2 ) = m 2 of degree 2. The maximum size of a clique for a given partial geometry is not known, 
but, depending on the values of r, 7 and k, they might form a significant trapping set since, 
from the above argument, assuming a clique of size k exists, t j k = p — k. A clique of size > 2 
is not possible for GQs since it implies forbidden triangles. 
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The cases above are extreme but serve as demonstrations that the bounds on the sizes of the 
trapping sets might be sharpened for the case of partial geometries and GQs. 

VIII. Expansion properties of the graphs of partial geometries 

The notion of expansion in coding theory originated in the work of Tanner 11451 (under the 
name of the strongly related concept of superconcentrators). That work also discovered the 
important relationship between the second largest eigenvalue of the adjacency graph of the code 
and its expansion properties, a relationship that has been widely exploited by mathematicians 
and computer science researchers since that time. These issues are examined further here. 

Codes that can be encoded and decoded in linear time can be constructed by means of codes 
with suitable expansion properties ll39l . the expander codes, and such codes and their graphs 
have received considerable attention over the past two decades. The codes were not LDPC codes 
and their decoding algorithm was not a message-passing one. However, it can be argued [f45l 
that expansion properties are also of interest in belief propagation decoding algorithms for LDPC 
codes since the notion of graph expansion can be interpreted as a measure of connectivity and 
randomness of the graph, desirable properties for efficient decoding properties. Tanner used such 
arguments in his seminal paper p4|. In ll46l he showed that the minimum distance properties of 
graph-based codes with small second eignevalue relative to the largest eigenvalue were also good. 
However, it is in the work of Burshtein and Miller Cl that a more direct relationship between 
expander graphs and the performance of message-passing decoding algorithms is established. 
That work uses expander-based arguments to establish that for sufficiently long block lengths, 
once a message-passing algorithm corrects a sufficiently large fraction of errors, it will eventually 
correct all errors. The argument considers Gallager hard and soft decoding algorithms but will 
be applicable for a wider class as well. 

Since the second eigenvalue property is important for codes from several points of view, it 
seems worthwhile to examine the expansion properties of the graphs from partial geometries, a 
problem that is addressed in this section, with the thought of using this property as a possible 
distinguisher of codes for further examination. The extent to which these arguments are viable 
would have to be confirmed with simulation. 

The eigenvalues of graphs associated with partial geometries are well known and the con¬ 
tribution of this section is to modify the known results to the purpose of interest, namely the 
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expansion properties of graphs from these geometries. A brief overview of the key ideas involved 
is given. 

The majority of works on this topic consider regular graphs, graphs with each vertex having 
the same degree. Our interest is solely in the biregular bipartite case, defined as bipartite graphs 
with n left vertices, each of constant degree c and m right vertices, each of degree d. When m 
and n are understood, we refer to this as a (c, d) -biregular bipartite graph. The recent work of 
[fT 8 l is an important contribution to this problem as it is the first work known to the authors 
that considers the expansion properties of these graphs, apart from the original seminal work 
of Tanner |[45l . Large classes of such graphs have been considered in the literature as Tanner 
graphs associated with partial geometries and other combinatorial configurations. Much of the 
work on LDPC codes is in fact on such biregular graphs. 

The notion of graph expansion is defined in several related ways. The following definition 
m will be used since other results of that work will be of interest. The same notion was 
used in the original work of Tanner 11451 . which also considered biregular bipartite graphs. Let 
a general graph Q = (V, E) (not necessarily biregular bipartite) with set of vertices V and 
edges E be a connected graph with no self loops. Define the boundary <5X of a subset X C V 
with | X |< a | V | as the set of neighbors of X. Define the expansion coefficient c(a ) for 
| X |< a | V |, for some positive fraction a, of the graph by: 

^ • I *X | 

c(a) = mm . r— 

0 ^xcv mm{| X |, | V\X |} 

for set X C V. Interest is often in the case where a < 1/2 where the minimization in the 
denominator is not necessary. 

In the case that Q = (V, E) is a fc-regular (each vertex is of degree k), | V |= n), and 
let A = \(iij]o<i,j <n be its n x n adjacency matrix (i.e., its point-point adjacency matrix where 
dij = 1 if vertex v t is connected to vertex vf). The matrix A will have n real eigenvalues which 
are listed in decreasing order 

do ^ AL ^ ' A 


and n orthogonal eigenvectors. 

If Q is k -regular then /i 0 = k and p n _ 1 = —p 0 if and only if it is bipartite. It can be shown 
that for any family of /.'-regular connected graphs with number of vertices tending to infinity 



29 


will have 

lim inf /i\ > 2 \Jk — 1 . 

n=^oo 

A finite connected fc-regular graph will be called Ramanujan if all of its eigenvalues other than 
±k satisfy the bound 

H < 2 Vk — 1 . 

Such a graph is described as having a “small” second eigenvalue. The amount of literature on 
the search for graphs with this property is very large. 

As noted, the interest of much of the theory on LDPC codes is focussed on the construction 
of (c, d) -biregular bipartite graphs for use as Tanner graphs of the code and the recent work ltT 8 l 
is precisely on the expansion properties of such graphs. This section considers that work in the 
light of the graphs from partial geometries, all of whose eigenvalues are known. 

Let Vi be the set of c-regular VNs and V 2 the set of d regular CNs. Interest is in the case 
where the expansion in the previous definitions is for sets of VNs and X C Vi and its boundary 
is a subset of the CNs V 2 and bounds on the resulting expansion coefficients is in terms of the 
eigenvalues of the adjacency matrix. This case was included in the interesting recent work [U 8 l 
and several results from that work will be of interest here. That work defined the (c, d) -biregular 
bipartite graph Q to be Ramanujan if 

Hi{G) < Vc — 1 + Vd- 1 , 

a natural extension of the regular graph case by setting c = d, although not theoretically justified. 

For coding applications, the parity check matrix H is often the m x n matrix derived from a 
combinatorial configuration or finite geometry, where the m rows of H are identified with the 
blocks of the configuration or lines of the geometry and the columns of H with the points. It 
can be described as the line-point adjacency matrix of the structure. The work of llT 8 ll refers to 
this matrix as the transfer matrix. 

In terms of this matrix H, the adjacency matrix of the configuration or geometry which is 
n x n, a point-point adjacency matrix, (and not using the blocks as graph vertices) is easily seen 
to be 

A 1 = H T H-7l, (11) 

where there are 7 lines intersecting a point and I is a n x n identity matrix. 
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The related adjacency matrix A which includes both the left vertices of the Tanner graph 
(VNs of the code) and the right vertices (CNs of the code), is of the form 


A = 


O H T 
H O 


( 12 ) 


This is the matrix and graph whose expansion properties are of interest. Of course in the case 
of a (c, d) -biregular bipartite graph, the rows of the mxn matrix H have weight d and columns 
have weight c. Note that 


A 2 


h t h o 
o hh t 


(13) 


As discussed in [0, fl2Dl . the eigenvalues of the adjacency matrix of graphs of partial geome¬ 
tries are well known. Such combinatorial structures correspond to two class association schemes 
and strongly regular graphs and the eigenvalues of the adjacency matrix of the graphs of the 
PaG( 7 ,p, 8) Ai of (fill) are given by: 


(p — 1)7 multiplicity 1, 

p-l-S multiplicity (14) 

-7 multiplicity (p ~ 1)( ^ ) | ( 7_~i 1 ^)~ 1)+a) • 

It is clear |[2TI that the eigenvalues of H T H are those of the matrix A, of (fill) , with 7 
added to each eigenvalue, and with the same multiplicities, i.e., the eigenvalues of H 7 H are 7 p, 
7 + p — 1 — 5 and 0 . 

It is noted in j45l that if x is an eigenvector of H 7 H with eigenvalue A 7 ^ 0 then x' = Hx 
is an eigenvector of HH 1 with the same eigenvalue. Thus H T H and HH 1 have the same 
eigenvalues (with different multiplicities). It is not difficult to show that if the eigenvalues of 
the matrix of interest A in (fl2l) are {p^ i = 0,1,..., n — 1}, then the eigenvalues of H r H 
are {/r 2 , i = 0 ,1 ,..., n — 1} and the expansion properties of the graph are determined by the 
“second” eigenvalue 

As A is a real symmetric matrix we have the rank of A, 774 , is equal to the rank of A 2 
and from (fl2l ) . ta = 2rn- The eigenvalues of A 2 are those of H T H which are the same (in 
magnitude, not necessarily in multiplicity) as those of HH 3 . Thus, arguing in the reverse from 
the previous paragraph, the eigenvalues of the matrix A of (fill) are ±^fp, ±^7 + p — 1 — 5 
and 0. Thus the second eigenvalue is y/'y + p — 1 — 5. 
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The expansion coefficient of any (c, d )-biregular bipartite graph is given by lfT8l . Ii45l 


c(a) > 


acd + n1( 1 — a) 


( 15 ) 


and the importance of the magnitude of the second eigenvalue is seen, i.e., the smaller the size 
of Hi the larger the expansion of the graph. 

Such information is used in ||46l and lt2Hl to determine bounds on the minimum distance of 
the codes with parity-check matrix H. The interest here is in the expansion properties of the code 
with parity-check matrix H. Since the second eigenvalue of the matrix A is V 7 + p — 1 — 8, 
this suggests that the graphs of partial geometries with 8 as large as possible, while still being 
partial geometries, will have the best expansion properties, i.e., 5 = 7 — 1 corresponding to nets. 
Notice that the two classes of partial geometries constructed in this paper are two classes of nets. 
For the first class, 7 = p = P and 8 — p — 1, and for the second class, 7 = p = t and 8 — t — 1. 
Hence, they have good expansion properties. While the parameters of the partial geometry, 7, p 
and 8, are not independent, the conclusion is interesting and worthy of further consideration. 

For interest, to conclude this section, the ratio of is compared for biregular bipartite 

graphs and k- regular graphs. From the discussion, for /.'-regular graphs, for the Ramanujan case, 


Mi/mLx = 2 (k - 1 )/k 2 « 2/k 


while for (7, p )-biregular bipartite graphs from partial geometries the ratio is 

,,2/,,2 _7 + p — (5 + l)_l , 1 (£+1) 

Ml/ Mmax ^ ' 

ip p 1 ip 

If parameters can be chosen so that p « 7 & k and 8 small, the two ratios are similar. For larger 
values of 8, however, such as for nets, the ratio of the eigenvalues of the partial geometry is 
approximately 1/p, suggesting graphs with large p and 8 would have better expansion. 

IX. Conclusion and Remarks 

Several aspects of codes derived from partial geometries have been considered. New results 
on codes from partial geometries were given, including an interesting characterization of them 
in terms of arrays of cyclic permutation matrices and two new and simple constructions. The 
trapping sets of codes from partial geometries were also investigated using the geometric prop¬ 
erties of their constructions. Finally, comments were given on the expansion properties of graphs 
from partial geometries. 
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